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We present the first ab initio calculation of a radiative transition of a hadronic resonance within 
quantum chromodynamics (QCD). We compute the amplitude for tttt —>■ vrq*, as a function of the 
energy of the tttt pair and the virtuality of the photon, in the kinematic regime where tttt couples 
strongly to the unstable p resonance. This exploratory calculation is performed using a lattice 
discretization of QCD with quark masses corresponding to m.,r ~ 400 MeV. We obtain a description 
of the energy dependence of the transition amplitude, constrained at 48 kinematic points, that we 
can analytically continue to the p pole and identify from its residue the p —>■ vrq* form factor. 


a. Introduction: The electromagnetic transitions 

of the nucleon into unstable resonant N* excitations are 
primary tools in the experimental study of nucleon struc¬ 
ture and spectroscopy [1]. These processes give us in¬ 
sight into the mechanisms that lead to the formation of 
the low-lying and excited hadrons from the basic quark 
and gluon building blocks of quantum chromodynamics 
(QCD). It is crucial to have a complementary theoretical 
program that connects physically observed transitions to 
QCD. One major challenge in studying these transitions 
is their resonant nature, where the N* excitation decays 
rapidly to asymptotic scattering states composed of two 
or more stable hadrons. To investigate these processes 
within QCD, one needs a nonperturbative framework 
that can accommodate resonant behavior, and presently, 
lattice QCD is the only available tool to evaluate such ob¬ 
servables while making only controlled approximations. 
Its implementation for the determination of properties of 
hadron resonances is still at an exploratory stage, and 
in this work we will extend the exploration into a new 
area with the first calculation of a radiative production 
amplitude of an unstable hadronic resonance from QCD. 

Before attempting the more complicated baryonic case 
of 7 * TV —)■ N* —>■ Ntt, we will consider a simpler problem 
featuring only mesons, Try* — > p — >■ tttt, which in addition 
to serving as the first of a new class of observables to be 
studied, is itself of significant phenomenological interest. 
The amplitude for this process is related to the hadronic 
contribution to the anomalous magnetic moment of the 
muon [2, 3], the chiral anomaly [4, 5], and the p -T Try ra¬ 
diative decay rate [6, 7], and appears in meson-exchange 
models of nuclear structure [8]. The p -T Try* transition 
has been previously studied using lattice methods (see, 
for example. Refs. [9-11]), but prior to this work the p 
has always been treated as a stable hadron, incapable of 
decay to tttt, in contrast to how it appears in experiment. 
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This approximation, which is uncontrolled for light quark 
masses such as those used in Refs. [10, 11] is removed in 
the present work. 

The perturbative nature of quantum electrodynamics 
ensures that to an excellent approximation the tttt -t Try* 
amplitude can be obtained from matrix elements of the 
electromagnetic current, 

= (out; Tr,P^|j//o|in; TTTT, £ = l), (1) 

where the tttt state has been projected onto an £ = 1 
partial wave, and where and Pn- are the 4-momenta 
of the TTTT and tt states, respectively. We will determine 
this amplitude as a function of the c.m. frame energy 
of the pion pair E*^ and the virtuality of the photon, 
= —{Ptt — P-kttY, by evaluating correlation functions 
using lattice QCD. 

Lattice QCD calculations are performed in a finite, dis¬ 
cretized Euclidean spacetime, and this introduces three 
length scales into the theory: the lattice spacing, a, and 
the spatial (P) and temporal (T) extents of the volume. 
Eor studies of stable hadrons not featuring heavy quarks, 
provided rriT^T, ^ 1 and a ^ 1 fm, the typical 

length scale associated with hadrons, these approxima¬ 
tions introduce only small and controllable systematic 
errors. 

The restriction to a finite volume in space prohibits the 
definition of asymptotic states, making the relationship 
between few-body observables obtained via lattice QCD 
and the scattering amplitudes of infinite-volume QCD 
somewhat nontrivial. As has been extensively explored in 
the literature, scattering amplitudes of two-body [12-19] 
and three-body systems [20-23] can be determined from 
the spectrum of eigenstates of QCD in a finite volume. 
Such spectra can be obtained from two-point correlation 
functions within lattice QCD, and the energydependence 
of hadron scattering amplitudes can be inferred - by an¬ 
alytically continuing these amplitudes to complex values 
of the scattering energy, poles can be found, with the 
pole positions providing the mass and width of hadronic 
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resonances. For an example see the recent determination 
of kaon resonant excitations in coupled-channel nK, rjK 
scattering [24, 25], and the p resonance for lighter quark 
masses where TnT,KK can be coupled [26]. 

The extension of the formalism to the case where an 
external (e.g. electroweak) current causes a transition 
from a single stable hadron to a pair of hadrons was 
presented by Lellouch and Liischer. They demonstrated 
that one can constrain such an amplitude using hadronic 
matrix elements of the currents evaluated in a finite vol¬ 
ume [27]. Their work focused on the implications of this 
formalism for K —>■ tttt decays, where tttt is in an S wave 
(see Refs. [28-32] for numerical implementations), and 
this has been subsequently extended to other systems of 
interest [15, 16, 18, 33-37]. 

Recently, these ideas were extended to accommodate 
more generic processes featuring an external current [38, 
39], and in this work it was shown that the transition 
amplitude, , can be obtained from finite volume 

matrix elements of the vector current. 


\n 
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where TZ is the residue of the finite-volume two-hadron 
propagator, which depends on the tttt 4-momentum, the 
cubic irreducible representation (A,n.,r)i the lattice volume 
{L X L X L), and the tttt elastic scattering amplitude. 
The hadronic finite-volume eigenstates carry labels. A, 
which indicate in which irreducible representation, or “ir- 
rep”, of the reduced rotational symmetry of the cubic 
lattice they lie. We point the reader to Ref. [39] for a 
detailed derivation and definition of TZ. Equation (2) is 
an approximation of the result presented in Ref. [39] - 
we have ignored contributions due to mixing with higher 
partial waves, which are both kinematically and dynam¬ 
ically suppressed in the energy regime of interest. ^ 

b. Lattice QCD calculation: We use an anisotropic 
Symanzik improved gauge and Clover fermion actions 
with Nf = 2 -|- 1 dynamical fermions. The quark masses 
are chosen so that m,, ~ 400 MeV [41], and we use a 
spacetime volume of {L/qsY x {T/at) = 20^ x 128, where 
Os and at are the spatial and temporal lattice spacings 
with Qs/at = 3.444(6) and Og « 0.12 fm. In Ref. [42] 
it was demonstrated that exponential corrections associ¬ 
ated with the finite volume of this lattice lead to subper¬ 
cent corrections. 

We construct three-point correlation functions using 
the technology presented in Ref. [9]. We use variation- 
ally optimized tt and isospin-1 tttt operators, and 


^ In Ref. [40] it was demonstrated that the £> 3 mr scattering 
phase shifts are consistent with zero. This assures one that the 
only resonance present that couples to the 1 = 1 tttt channel in 
this kinematic regime is the p-meson and so the £ = 1 transition 
amplitude is the dominant contribution. 


respectively, that have been subduced to the de¬ 
sired irrep. A, of the appropriate little group of the octa¬ 
hedral group [43]. These operators have been previously 
obtained in the determination of the spectrum from two- 
point correlation functions [40] . Inserting the vector cur¬ 
rent we have three-point functions, 

(0|fli^-l(At, P^):7'^(t,P^-P^^)4V^(0,P,^)|0), 

( 3 ) 

and we will present results extracted from correlation 
functions computed with Euclidean time separation. 
At = 32at, excluding Wick contractions where the cur¬ 
rent couples to a disconnected quark loop is the 

tree level improved Euclidean vector current [9], which 
is renormalized by insisting the tt form factor be 1 at 
= 0 GeV, giving a multiplicative renormalization of 
Zy = 0.833(9). By inserting a complete set of finite 
volume QCD eigenstates in Eq. (3), and evolving the 
operators to the origin of Euclidean time, one can de¬ 
termine P^., -P-n-Tr; A,r 7 r)^ from the time 

dependence of the correlation function [9]. 

Consideration of various momenta, P = ‘^\nx, ny, Uz], 
allowed by the periodic boundary conditions, leads to 
determination of the matrix element at 48 distinct kine¬ 
matic points. Eight different discrete E*^ values feature, 
corresponding to the finite-volume eigenstates of tttt in 
various irreps, and discrete values of photon virtuality in 
the range —0.4 < (Q/GeV)^ < 1 are sampled. 


c. The transition amplitude and p —?> Try* /orm fac¬ 
tor: The TTTT —>■ TTTT P-wavB elastic scattering amplitude, 
expressed via a phase shift, 6i{E*^), was determined 
from the lattice QCD finite-volume spectrum in Ref. [40] . 
With this in hand we may evaluate TZ in Eq. (2) and de¬ 
termine the transition amplitude from the finite-volume 
matrix elements. 

The transition amplitude is a Lorentz vector, and it 
has a kinematic decomposition, 

n/U _ p p cr,(\ PI 2 A 

^7777 777* — ^ -^77,I/-‘7777,aCy3 A7r77 1 ^ 7777 J ^ ^7777,777* 5 

( 4 ) 


where Q"^) is a Lorentz scalar and ep is the 

polarization vector of the P-wave tttt state with A,r,r being 
its helicity. 

In the energy region we consider, the transition ampli¬ 
tude will be sharply peaked due to the p pole. Defining 
a p —>■ Try* form factor, Ft^p{E*^,Q'^), we may write the 
amplitude, 


•^irTT.'K'y* Q^) 


= p.p(p;,,q2) 


8Tr 






( 5 ) 


^ These exactly vanish in the SU{3) flavor limit and are expected 
to be suppressed for the quark masses used. 
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FIG. 1. The points appearing in the upper panel depict the 
form factor determined from lattice QCD for three tttt energy 
levels. The index n labels the order in which the state appears 
in the spectrum. Also shown are the fits of the form factor 
performed using Eq. (6) and evaluated at the three tttt ener¬ 
gies. The corresponding P-wave phase shift and tttt energy is 
shown on the lower panel. 




which is proportional to the elastic tttt — >■ tttt scat¬ 
tering amplitude sin^i with = 

In Refs. [34, 38] it was demonstrated 
that this parametrization is consistent with the analytic- 
ity and unitarity constraints required in scattering the¬ 
ory. The presence of the energy-dependent p —>■ tttt 
strong decay width, Ti, can be understood in the con¬ 
text of an effective field theory where the rescattering 
of the final tttt states is mediated by a fully dressed p 
resonance (see the appendices of Ref. [38]). 

The derivative of the phase shift, , appears in IZ 
- to compute it we use a sensible parametrization for 
the relativistic Breit-Wigner function [40]. 

In Fig. 1 we present the computed form factor for three 
of the eight tttt energies studied, using the Breit-Wigner 
parametrization for the phase shift, where we observe 
that both spacelike and timelike kinematics are sam¬ 
pled. It is evident that F.,^p[E*^,Q'^) has only a mild 
dependence on E*^, with the sharply peaked resonant 
behavior having been captured by the sin(5i(P*^) factor 
in Eq. (5). 

To analytically describe the E*^ and dependence 
of the form factor we introduce an ansatz, 

l + a.<3- + A(£g-m§) + 

-b as exp [-aeQ^ - P 2 (.Ell - ml)] 

+ P,iE:l-ml)+^,iE:t-mt), ( 6 ) 

where the parameters a^ and /3i are to be fitted and the 
constant mg is fixed to 2.1805 to coincide with the real 
part of the p mass. To fit the form factor, we vary the 
form being used by setting a subset of these coefficients to 
zero and thus consider over 15 different fit functions. We 
also consider fits where the points in the timelike re¬ 



FIG. 2. The upper panel shows the real part of the form 
factor determined in this work [orange band] evaluated at the 
p pole, Eplm,r = 2.1805(32) — 0.0151(5)i. This is compared 
with the value obtained in Ref. [9] [green band], where the p 
resonance is QGD stable, and the experimentally determined 
pn photocoupling [6, 7]. The lower panel shows the previously 
undetermined imaginary component of the form factor. 


gion are excluded. From all fits performed, we retain only 
those that have a y^/DOF < 1.5, and we find that no fit 
lacking i?*,r"dependent terms satisfies this. The bands 
shown in Figure 1 reflect the parametrization variation 
as will the uncertainties on all quantities quoted below. 

With an analytic description of the E*^ dependence 
of form factor, we may analytically continue to the p 
pole at E^^ = [2.1805(32) - f 0.0151(5)] The 
dependence of the resulting form factor is shown in Fig. 2, 
with the small imaginary part reflecting the fact that the 
p resonance in this calculation is unstable, but with a 
small hadronic width - as the pion mass is decreased in 
future lattice calculations [26], the width will increase 
and with it the imaginary part of the form factor. 

The transition amplitude, 47r7r,7r7*) follows from 
Eq. (5,) where the phase is fixed up to an overall sign by 
Watson’s theorem to be the tttt —>■ tttt phase shift. The 
remaining sign only has meaning in comparison to other 
transition amplitudes, and consequently, we need only 
present the absolute value of In Eig. 3 we plot 

w,n.|^ 7 r 7 r.ir 7 * I as a function of E\^ for two values of . 
This figure illustrates that as the tttt energy approaches 
the p pole, the transition amplitude is dynamically en¬ 
hanced by the resonance as one would expect. The res¬ 
onant behavior, as a function of arises solely from 
the IZ factor in Eq. (2); it is not due the parametrization 
in Eq. (5) which simply serves as the definition of the 
form factor 

From Eq. (4), one may readily obtain the Tr+y —>■ 
cross section in terms of the reduced amplitude, 
evaluated at = 0, 


(T(7r'^7 —>■ = 


^TTTT Qnj 


A 


•TTTT,777' 


iKl,0) , ( 7 ) 


where q*^ are the c.m. relative momenta. In Fig. 4 
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FIG. 3. The top panel shows m.,r|.47r7r,7r7* | as a function 
of the c.m. tttt energy. This is determined for two different 
values of Q^/GeV^ = 0,0.803. For comparison, in the lower 
panel we show the absolute values of the elastic £ = 1 tttt 
amplitude, 


FIG. 4. The top panel shows the vr+y —>■ tt+tt® cross section 
as a function of the tttt c.m. energy. The lower panel shows 
the elastic £ = 1 scattering cross section. One observes near 
the resonance the enhancement of the 7r^7 —>■ tt+tt® cross 
section. 


we plot this as a function of the c.m. energy. That 
the peak cross section for ^ 400 MeV is significantly 
larger than phenomenological parametrizations of the 
physical cross section [44, 45] can be easily understood: 
near the resonance we have 


lim —>■ oc 

E*^^mp 


< 1 % 0 ) 

ml Fi (mp) 


and the F^p(rnp,0)/ml ratio we find to be approx¬ 
imately 60% of the experimental value, and we expect 
this to vary only slowly with changing quark mass. On 
the other hand, the width of the p resonance when 
400MeV, 12.4(6) MeV [42], is approximately 12 
times smaller than the experimental width [46], scaling 
as expected for an approximately quark mass indepen¬ 
dent coupling, gpTT-K, with reduced P-wave phase-space. 
This suggests that as future calculations are performed 
with quark masses closer to their physical values, and as 
the p resonance becomes broader [26], the Tr+y —7r+7r° 
cross section will decrease by an order of magnitude. For 
comparison, in Fig. 4 we plot the £ = 1 elastic 

cross section, whose factor of 5 kinematic enhancement 
with respect to the experimental determination (see for 
example. Ref. [47]) can be understood by the 1/g*^ de¬ 
pendence in the vicinity of the resonance. 

d. Final remarks: We have presented the first de¬ 
termination of a resonant radiative transition amplitude 
from QCD. This exploratory study of tttt — )■ Try*, al¬ 
though performed with unphysically heavy light quarks, 
serves as a proof of principle that hadronic transition pro¬ 
cesses involving resonating few-body states can be rigor¬ 
ously studied using lattice QCD. We have demonstrated 
how from this amplitude, by analytically continuing to 
a pole in the complex energy plane, one may obtain the 
p —>■ Try* form factor where the p is treated as an unsta¬ 
ble resonance, and have also obtained the Tr+y —>■ tt+tt^ 


cross section, and discussed how we expect the results to 
change in future calculations using lighter quark masses. 

Closely related techniques can be implemented in fu¬ 
ture studies of hadron structure and weak decays. As 
well as the obvious extension into the baryon sector, 
y*V —>• N* Ntt, there are processes important 
for testing the limits of the standard model such as 
B —)■ Ktt£'^£~ [48, 49], where the Ktt system is known 
to resonate. 

Having demonstrated in this work the feasibility of 
studying radiative transition of two-body hadronic res¬ 
onances directly from QCD, future studies will focus on 
the extension of this work. The technology for studying 
transition amplitudes with any number of open two-body 
states has been already developed [38, 39] and here we 
have tested it in the case where there is only one chan¬ 
nel open. Future calculations will accommodate simi¬ 
lar processes involving resonances that decay strongly to 
more than one hadronic state, for example iFy* —>■ K* —>■ 
Kir/Kr] [24, 25] and Try* — p* — >■ tttt/KK [26]. Further¬ 
more, given the recent and exciting theoretical develop¬ 
ment for the study of three strongly interacting parti¬ 
cles via lattice QCD [20-23], we can also expect electro¬ 
magnetic transition amplitudes involving three or more 
hadrons (e.g. Vy* —)• N* Ntttt). ^ 
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